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wave coupling between the local moments and the conduction electrons on the ladder. By changing 
the relative strength of the two Kondo interactions we examine the evolution of the system from a 
conventional Kondo insulator with a singlet at each site to a new kind of semimetallic state formed 
by overlapping of Zhang-Rice-like singlets. The DMRG is used to study how the spin and charge 
correlation functions evolve between these two regimes. 
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I. INTRODUCTION 



The interaction between localized spins and delocalized 
conduction electrons plays a central role in the physics of 
heavy fermion materialsjj Kondo insulatorstm and many 
other highly correlated electron systems. The Kondo lat- 
tice model (KLM) is probably the simplest theoretical 
model that describes this physics. 

Considerable insight into local moment materials has 
been obtained by studying the one-dimensional Kondo 
lattice.Q Although some properties are intrinsic to-pne di- 
mensional models, such as the nesting instabilities^, other 
aspects such as the insulating behavior at half filling, and 
the Fermi surface expansion by the local momenta!! are 
also inherent to three-dimensional systems. Recently, 
the theoretical study of one-dimensional models has also 
proved a most useful tool in the analysis of an increasing 
number of low-dimensional materials. The comparison of 
experimental findings with exact numerical results helps 
to evaluate the adequacy of basic models and to quantify 
the relevance of additional effects such as the importance 
of the orbital electronic degrees of freedomM 

The earliest claims of quasi-one-dimensional behav- 
ior in Kondo lattice, intermetallic compounds date back 
to the seventies ,BEJ but only recently, with the syn- 
thesis of the quasi-one-dimensional molecular conductor 
(DMET)2FeBr4, has a clear-cut experimental realiza- 
tion of a ID Kondo lattice been foundO This compound 
displays features of a strongly correlated electron system. 
It shows a metal- insulator transition at 200 K, which 
appears to be driven by the interaction between the 7r- 
conduction electrons f>£ the DMET donor sheets and the 
localized Fe 3+ spinsJi3 

A class of heavy fermion systems, known as 
Kondo insulators, display insulating properties at low 
temperatureslm One of the interesting aspects of a 
Kondo lattice is the emergence of insulating ground- 



states. The ground state of the conventional Kondo lat- 
tice model is a non-magnetic insulator at half-filling, with 
a gap in both spin and charge excitation spectrum. □ In 
this way, the half-filled Kondo lattice provides a simple 
model for the Kondo insulator, where the coherent Kondo 
effect forms a highly renormalized band-insulator. How- 
ever, the smallest gap Kondo insulators, CeNiSn and 
CeRhSb, do not fit into this simple scheme: they appear 
to develop gapless excitations.E 3 ] The narrow-gap Kondo 
insulators can be regarded as an example of a half-filled 
Kondo system with a gapless ground state due to the 
presence of interactions that-are absent in the conven- 
tional Kondo lattice model.tfJIij 

A simple model which displays a metallic ground state 
at half-filling is the "d-wave Kondo model". This is a 
variant of the conventional Kondo model where each lo- 
cal spin couples through the exchange interaction with a 
d-wave orbital formed by the nearest-neighbor electrons. 
This interaction favors the formation of Zhang-Rice-like 
singletalj between each spin and the nearest conduction 
electron orbitals. The superposition of these extended 
singlets may establish a new kind of semimetallic state, 
which may be relevant to explaining the metallic phase of 
{D M ET)2F eBr^. Another interesting feature of the d- 
wave Kondo model is the introduction of spin-dependent 
hopping processes, where an electron is able to jump to 
the next-nearest-neighbor site by flipping the spin of the 
local moment. This kind of interaction plays a role in 
other organic compounds_-such as the magnetic graphite 
intercalation compounds.ED 

Motivated by the need to understand how a metal- 
insulator transition arises from pure exchange interac- 
tions between the itinerant and localized electrons, this 
paper examines a Kondo model with two spin-exchange 
channels sharing a single conduction band: an on-site 
s-wave (Fig. fy) and a nearest neighbor d-wave cou- 
pling (Fig. |^). Electron-electron interactions, such as 
Hund's interactions, can induce the opening of new spin- 
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exchange channels between the local moments and the 
conduction electrons. c3~E3 For an impurity magnetic ion, 
the Kondo effect develops exclusively in the strongest 
screening channel due to the local symmetry which pre- 
serves the channel quantum number of scattered elec- 
trons. In a lattice, however, a conduction electron can 
change symmetry channels as it moves from one spin site 
to another. Therefore, the development of a Kondo ef- 
fect in one channel no longer excludes the possibility of 
a Kondo effect developing coherently in other channels. 

Our goal is to determine whether there is a range of 
exchange coupling for which such a two-channel Kondo 
model displays gapless excitations at half-filling. If this 
is the case, the competition between the d- and s-wave 
Kondo couplings might induce new intermediate phases 
between the d-wave Kondo ground state (Zhang- Rice sin- 
glets at each site) and the conventional Kondo ground 
state (s-wave Kondo singlets). It has been argued that 
even a composite superconducting phase could appear 
in the KLM if the local moments interact with a single 
conduction band via two orthogonal scattering channels 
with the same spatial parity.EU 

In highly correlated electron systems, it is essential 
to go beyond the mean-field level to understand their 
low-energy properties. We apply the .density matrix 
renormalization-group (DMRG) methodca to calculate 
the energy spectrum and ground state correlations of this 
two-channel KLM. This real-space technique has proved 
to be remarkably accurate for one-dimensional systems 
such as the Kondo and Anderson lattices£j~Ej To check 
the character and energy range of the elementary excita- 
tions we have also performed a variational calculation in 
the d-wave limit. 

Section II introduces the Hamiltonian and the param- 
eters of our model, the scaling law used to analyze the 
results and the numerical technicalities of the calcula- 
tion. We present our results for the ground state energy 
and the spin and charge excitation gaps in section III. In 
Section IV, we discuss the spin and charge correlations of 
the ground state and some characteristics of the excited 
states. Section V is dedicated to a comparison between 
our results and a variational calculation and finally, sec- 
tion VI presents our conclusions. 



II. MODEL 

The minimal one dimensional model with two Kondo 
exchange channels of different symmetry and the same 
parity is a ladder model. The exchange interactions be- 
tween the local moments and the electrons are an on-site 
s-wave (Fig. [I]) and a nearest neighbor d-wave coupling; 
(Fig. ||). Unlike other multichannel Kondo models 
our model involves a single conduction electron band 
and Kondo singlets can form for any value of the cou- 
plings. Electrons and local moments are located at sites 
= ix + ky, where k = 0, 1. We impose open bound- 



ary conditions in the x direction. Along the y axis, the 
system can be considered a periodic lattice with a period 
of two lattice spacings. The Hamiltonian is: 

H = H + H s +H d , (1) 

where Hq is the electronic kinetic energy given by, 

H = -*5>*I la V& + V>L +1 VL + h-c), (2) 

where t denotes the hopping matrix element, and the op- 
erator ipil creates a conduction electron on site R% with 
spin a. H s is the conventional s-wave Kondo interaction: 

H s = J s J2{%-Tp£^^} (s- channel), (3) 

i,k 

where J s is the strength of the antiferromagnetic s-wave 

coupling, and sf = ipil(crap/2)ipip denote the spin of 
the local moment and the conduction electron at position 
R%, respectively. 




FIG. 1. Ladder illustrating the electronic hopping and the 
s-wave Kondo interaction between each local spin (arrows) 
and the conduction electron at the same site (spheres). 

The last term, H dl represents the d-wave Kondo chan- 
nel: 

H d = J d J2{S--d^4 p } (d- channel), (4) 

i,k 

, ,kf (^ k +la+VA k -l a ~7V'i ( a + 1)t ) 

where = Af==f ■ ( 5 ) 

The operator d*^ creates an electron with "d-wave sym- 
metry" in the "7" Wannier state at site Rf. For 7 = 1 
this state has pure d-wave symmetry but since the one- 
dimensional ladder breaks the 90° rotation symmetry of 
the lattice, we explore a more general class of states by al- 
lowing 7 to become a free parameter in the model (in this 
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way our "d-channel" involves an admixture of d and ex- 
tended s-wave character) . Note that the d-wave spin ex- 
change involves spin-dependent hopping processes where 
the electron hops between next-near-neighbor sites via 
the exchange of spin with the local moments. 




FIG. 2. Ladder illustrating the d-wave Kondo interaction 
between each local spin and the three near-neighboring elec- 
trons. Note that the coupling is Jd along the chains and —~/Jd 
along the rungs of the ladder. 



Our interest lies in the half-filled case where the total 
number of conduction electrons is equal to the number 
of lattice sites. Finite size studies on our model are fa- 
cilitated by using the strong-coupling limit to reduce the 
feasibility of antiferromagnetic instabilities, although, as 
we discuss later, we can not completely eliminate them. 
Thus, we use the following values of the parameters, 
t = 0.01 and |J| = y/ Jf + J'j = 1. A finite hopping 
t is required to ensure the smooth distribution of the 
electronic density. We have studied three values of the 
transverse exchange parameter, 7 = 2, 1, 0.5, and several 
values of the ratio between the d- and s-wave exchange 
couplings, J d /J s . |-, 

We use the finite-chain DMRG algorithmic up to lad- 
ders of 12 rungs, using up to 100 optimized states for 
each block. In contrast to the infinite system method, 
the finite system algorithm gives more accurate results, 
but an asymptotic extrapolation to the infinite-,system is 
necessary. We use the following scaling law,E3 which is 
appropriate for a massive mode, 



A(L) = yA2(L = co)+^__^ ( 6 ) 

where A(L) is the gap for a system of length L. This al- 
lows us to define a correlation length, £ = v/A(L = 00), 
for the spin and charge modes of the system. In the 
DMRG calculations, the truncation of the Hilbert space 
might lead to deviations from the above asymptotic 
behaviorO The deviations can be substantial for small 
gaps, especially when the dimension of the Hilbept space 
used in the computation is not sufficiently large.E3 

In our computation, the largest truncation errors are 
of the order of 10~ 7 for the conventional s-wave limit, 
while for the biggest Jd values considered, the errors are 



as big as 10~ 3 . As the ratio Jd/J s increases, the trun- 
cation error becomes larger since the Kondo screening 
cloud extends across substantially longer distances and 
the magnetic fluctuations grow. As a result, the accuracy 
of our calculation is significantly reduced in this limit. 



Ground state energy 




Arctg(J d /J s )/7i 

FIG. 3. Ground state energy per rung of the two-channel 
Kondo ladder model versus arctan ( Jd/ Js)/k for transverse 
exchange parameter 7 = 2 (dark gray symbols), 7=1 (black 
symbols) and 7 = 0.5 (light gray symbols). We have plotted 
the results for ladders of 5 (circles), 6 (squares), 7 (diamonds), 
8 (up-triangles), 9 (down-triangles) and 10 (stars) sizes. Note 
that the lines in this plot are only guides to the eyes. 



III. GROUND STATE ENERGY AND SPIN AND 
CHARGE EXCITATION GAPS 

Let us define the ground-state energy in a given spin-S 
subspace for a finite ladder with L- rungs by E g (L, N c , S), 
where N c is the number of conduction electrons. Given 
that we are interested in the insulating ground state, we 
will consider a half-filled band with one electron per site 
and a total of N c = 2L electrons. The spin gap is de- 
fined as the energy difference between the ground state 
and the lowest spin excitation that changes the total 
spin quantum number by one: A S (L) = E g {L, 2L,S — 
l)-E g (L,2L,S=0). 

Similarly, the charge gap is the energy difference be- 
tween the ground state and the lowest pure charge ex- 
citation, which changes the total carrier number by two 
and keeps the spin quantum numbers fixed, A C (L) — 
E g (L,2L + 2,0)-E g (L,2L,0). 
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Fig. H displays the energy of the ground state as a 
function of arctan(Jd/J s )/7r for the three values of 7 ex- 
plored in our calculation. Results for several lattice sizes 
have been plotted. The ground state energy increases 
with the strength of the d-wave coupling for any value of 
7 and reaches the highest values when 7 = 0.5. 

Let us now discuss the results for the ground state 
and excited state energies of the conventional Kondo lat- 
tice model, J 8 = 1 and Jd = 0. In the strong cou- 
pling limit, the energies can he calculated by using a 
perturbative expansion in tj J.a To second order in this 
small parameter, the ground state energy per rung is 
E g = -(3/2) J - (4i 2 /J) = -1.5004, the spin gap is 
A., = J - (20t 2 /J) = 0.998, and the charge gap is 
A c = (3/2)J-4t+{9t 2 /J) = 1.4609, where we have used 
the parameters of our model. Our numerical results are: 
E g = -1.50040, A s = 0.9981 and A c = 1.4614, where 
we have used the scaling discussed in the previous sec- 
tion. Given the almost complete agreement between the 
perturbation expansion and our numerical results, the 
accuracy of our computation in this limit is confirmed. 

In the d-wave limit, J s = 0, Jd — 1, the ground state 
energies for the three different values of 7 are the follow- 
ing: 



Charge and spin gaps 

y=2 



• 7 


— 2 E run9 ~ 

' ground 


-1.05 


• 7 


— 1 zp run 9 
' ground 


-0.93 


• 7 


-05 E rung 


~ -0.8 



These are reasonable values for the ground state energy 
in the d-wave strong coupling limit. For large 7, the 
system will tend to form singlets along the rungs of the 
ladder. The energy of such configuration is: E™ 1 ™ 



(~3/4)- For 7 = 2, the ground state energy is just 
Aground = — ^ ^he vame obtained in our calculation is 
rather close to this simple limit. In the opposite case 
of 7 = 0, let t = 0, there are four decoupled chains. 
The spins at each chain tend to form spin-singlets with 
their two neighboring electrons with exchange coupling 
Jd/2 = 1/2. Therefore, the ground state energy per rung 
is related with the ground state energy of four Heisen- 
berg chains: £™" 9 nd = 4(l/2)(l/4-ln2) = -0.886. The 
value obtained for 7 = 0.5 in our calculation is again close 
to this simple limit. The numerical ground state energy 
for 7 = 1 is between the upper bound of —1 and the lower 
bound of —0.886. Therefore, our result agrees with this 
simple physical argument. 

Figures ||, || and || display the charge and spin gaps 
versus arctan (J^/Js)/^ for 7 = 2, 1 and 0.5, respec- 
tively. Results for several lattice sizes have been plotted. 
The gap resulting from using the scaling law QG ) is also 
shown. In addition, the first column of tables j| and |H| 
displays the spin and charge correlation lengths, respec- 
tively, derived from the scaling parameters of equation 
(|U = «/A(L = oo). 
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FIG. 4. Charge gap (white symbols) and spin gap (black 
symbols) versus arctan (Jd/Js)/™ for a d-wave transverse ex- 
change parameter 7 = 2. We have plotted the results for lad- 
ders of 5 (circles), 6 (squares), 7 (diamonds), 8 (up-triangles), 
9 (down-triangles) and 10 (stars) sizes. Also, the charge (solid 
line) and spin (dashed line) gaps obtained using the scaling 
procedure are shown. 

The behavior of the spin and charge gaps for any value 
of the parameter 7 is similar. By increasing the d-wave 
coupling, the system evolves to a state with much smaller 
spin and charge gaps for any 7 = 2,1,0.5. Comparing fig- 
ures ||, H and [| it can also be seen that A s and A c in 
the d-wave limit increase as 7 decreases. So, a transverse 
exchange parameter of 7 — 2 is the most favorable value 
to develop a gapless region in the phase space. In fact, 
the smallest values of the gaps correspond to 7 = 2 and 
Jd/ J s = 10, where by using scaling (0) we are able to in- 



fer values of A s 



2 • 10~ 7 and A c 



0.034. Our 



belief is that an increase in values of the d-wave Kondo 
exchange coupling, which implies a derealization of the 
singlet formed between the spin and the conduction elec- 
trons, induces a reduction of the energy gaps. We think 
that for 7 < 1, the ladder sustains a spin and charge 
gap. However, for 7 > 1, the system develops gapless 
excitations at a finite value of the ratio Jd/J s - Unfor- 
tunately, the accuracy of the calculation is not enough 
to get precise results in the Jd/J s ^ 00 limit. We can 
not conclusively establish our expectation that a gapless 
region is present for any value of the ratio Jd/Js ft 10 
when 7 > 1. Also, we do not know whether spin and 
charge gaps drop to zero precisely at the same ratio of 
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Jd/Js or the spin excitation becomes gapless before the 
charge gap closes. 



Charge and spin gaps 
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Arctg(J d /J s )/7t 

FIG. 5. Charge gap (white symbols) and spin gap (black 
symbols) versus arctan (Jd/Js) /n for a transverse exchange 
parameter 7 = 1. Notations are the same as before. 



Other models also display gapless excitations at half- 
filling, such as the attractive Kondo-Hubbard model, 
where, in addition to the Kondo exchange, the conduc- 
tion electrons have an on-site attraction.Ea In the strong 
coupling limit of that model, there is a range of param- 
eters where both excitation gaps become equal and then 
drop to zero. However, in our model, although both 
excitation gaps are strongly reduced, the ratio between 
them grows when Jd/Js increases. This fact implies that 
the time scale of the collective spin fluctuations that de- 
stroy the antifcrromagnctic long-range order increases 
faster with Jd/Js than the characteristic time scale of 
the conduction-electron propagation. 



Charge and spin gaps 
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Arctg(J d /J s )/7t 

FIG. 6. Charge gap and spin gap versus arctan (Jd/ Js)/tt 
for a transverse exchange parameter 7 = 0.5. Notations are 
the same as before. 

Let us conclude this Section commenting on the nu- 
merical difficulties of the problem. The computation is 
much more involved in the d-wave limit (J s = 0, Jd = 1). 
The gap is, at least, two orders of magnitude smaller. 
Also, the correlation length in the ground state increases 
roughly by two orders of magnitude, and longer ladders 
are required to reach the convergent regime with respect 
to the chain length. The number of degrees of freedom in 
the Kondo lattice model is eight states per site. There- 
fore, the exact diagonalization can be carried out only 
for rather short ladders (five rungs in our case). 

In addition, as it is. the case for the one dimensional s- 
wave Kondo lattice^ we find that the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction increases in impor- 
tance as J s decreases. For small J s , large Jd, the system 
is so close to an instability that the numerical inaccuracy 
associated with larger ladder sizes is able to induce spuri- 
ous antiferromagnetic ordering of the lattice, because the 
states with and without long range order are very close 
in energy. 

Also, in the d-wave limit, the smallness of the hop- 
ping term makes difficult the redistribution of the local 
densities and, as consequence, the determination of the 
smooth and energetically lowest electronic states. Our 
approach to improve the numerical accuracy has been to 
make several self-consistent sweepings for each block to 
guarantee the smooth redistribution of the local densi- 
ties. 
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IV. SPIN AND CHARGE CORRELATIONS IN 
THE GROUND STATE. CHARACTERISTICS OF 
THE SPIN EXCITED STATES 

The DMRG method enables us to calculate the equal- 
time spin and charge correlation functions. We minimize 
the effects of the open boundary by evaluating correla- 
tions between the central site and other sites along the 
ladder. To clarify the nature of the correlations, we con- 
sider the pure d-wave exchange model with no hopping 
(t = 0). In this case, the ladder decouples in two discon- 
nected zig-zag chains. If only d-wave exchange is present, 
the spin sitting at position "0" in Fig. correlates only 
with the electrons at the black sites, not with the ones 
in the gray chain. Correlations between the central spin 
and electrons at the gray sites are induced by either s- 
wave exchange or hopping terms. We use black and gray 
colors to display the different correlation functions. For 
example, the spin-spin correlation (Sq ■ s*) between the 
central spin and the electron at site will be display in 
black if R\ belongs to the black chain and in gray in the 
opposite case. 
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sites decreases. When only d-wave exchange is present 
(J s = 0, Jd = 1) the correlation with the electrons in 
the black chain becomes dominant. The correlation with 
the electron along the rung reaches its highest value of 

(s*8**S> = -o.i6. 
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FIG. 7. The black sites couple with the spin through the 
d-wave Kondo channel. The gray sites couple through the 
conventional s-wave channel. 

In the calculation of the correlation functions, we 
have used the local spin symmetry of the ground state 
((pif) = (pij,)) and the fact that the ladder is half-filled. EHl 
We only discuss results for a d-wave transverse exchange 
value of 7 = 2, but the correlations for the other values 
of 7 exhibit similar characteristics. 

Figure |§| shows the spin-spin correlation between the 
local moment at the center of the ladder and electrons 
at different sites, (So • s*i), as a function of the electron 
position for several values of the ratio Jd/J s - It can be 
seen how the ground state evolves with increasing val- 
ues of Jd/J s from a conventional insulating state with 
strong on-site correlations to a new state with large spin- 
spin correlation along the rungs and longer correlation 
lengths. 

In the s-wave limit (top-left graph in Fig. ||) the 
value of the on-site correlation is as large as {S z qS z q) = 
—0.2499. The correlations with the nearest neighbors 
in this case are extremely small, (S z qS z i) ~ (S*qS x q) < 
0.001. This confirms the localized character of the ground 
state when only the on-site exchange is present. As Jd/Js 
grows, the correlation with the electrons at black sites 
gradually increases, and the correlation with the gray 
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Distance to the central site 
FIG. 8. Spin-spin correlation between the central spin and 
the electron at rung j as a function of the distance j. Corre- 
lations with electrons at gray (black) sites are shown in gray 
(black). The transverse exchange is 7 = 2. Values for sev- 
eral ladder sizes are represented with different symbols as in 
previous figures. 



We are also able to infer typical values of the spin cor- 
relation lengths by fitting the spin-spin correlation func- 
tions. An exponential fit, (S ■ Sj) — Cexp [— j/£ spin ], is 
appropriate for most of the ratios Jd/Js- However, for 
J d /Js > 10 a power law fit, (S ■ s}) = C/[(j /£ spin ) u + 1], 
more accurately describes the correlations along the gray 
chain, as it can be seen in the top graph of Figure (||). 
The transition from an exponential to a power law decay 
suggests the development of gapless spin excitations in 
the system, in agreement with our previous conclusions. 
Table | displays the correlation lengths between the cen- 
tral moment and the electrons at the gray (black) zig- 
zag chain, C^ay (£&fa™fc)> f° r different values of the ratio 
Jd/Js- The correlation length derived from the scaling 
of the excitation gap (£ s P m = v/A) is also shown. We 
notice that £ s ?> m and £|ray have similar values and both 

increase with the ratio Jd/Js- On the other hand, Cul™k 
is roughly constant for all the values studied. 
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correlation that becomes quite big in the d-wave limit, 
C(0) ~ 0.24. The correlation with the remaining sites 
in the ladder is smaller and negative. This large on- 
site charge-charge correlation is related to the non-local 
character of the singlet state in the d-wave channel, as 
we discuss in the next Section. 

We can also notice in Fig. [l(] that the correlation 
length for both zig-zag chains grows with Jd/J s - Table 
H shows the correlation lengths between the electron at 
the central site and electrons in the gray (black) zig-zag 

chain, £gray 9e (^wacfc 96 )' f° r different values of the ratio 
Jd/J s - Also, the correlation length derived from the scal- 
ing (£ char a e — v/A) is shown. All the correlation lengths 
grow with the ratio Jd/ J s , but their precise values are dif- 
ferent. As in the spin case, the charge correlation lengths 
have been obtained using an exponential fit in the region 
Jd/Js < 10 and a power fit for Jd/Js > 10. Middle and 
bottom graphs in Figure (j^) show a comparison between 
exponential and power law fits for Jd/Js = 10. 



1 2 3 4 5 6 

Distance to the central site 

FIG. 9. Comparison between exponential and power law 
fits to the DMRG results for the correlation as function of 
the distance to the central site for Jd/Js = 10 and a ladder 
of 10 sites. Top graph, middle graph and bottom graph dis- 
play the numerical data (circles), exponential (solid line) and 
power law (dashed line) fits for the spin correlation along the 
gray chain, the charge correlation along the black chain and 
the charge correlation along the gray chain, respectively 



TABLE I. Spin-spin correlation lengths for different values 
of the ratio Jd/Js (first column). The second column dis- 
plays the lengths derived from the scaling law. The third and 
fourth columns show the correlation lengths between the cen- 
tral local moment and electrons at the gray and black chains, 
respectively. 
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0.72 
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0.71 




0.88 
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3.17 


0.64 



Let us now discuss the characteristics of the charge 
correlation functions. Figure [To] displays the charge cor- 
relation, C(j) = (p(O)p(j)) - (p(0))(p(j)}, as function 
of j, the position of the second electron. The charge- 
charge correlation is practically zero everywhere for a 
pure s-wave exchange. This behavior agrees with the 
strong coupling limit of the Kondo model. The introduc- 
tion of a finite d-wave coupling induces an on-site charge 



a. 

V 

/\ 
o 

V 

I 

/\ 
3 

Cl 



Q. 

V 




0123456 123456 123456 

Distance to the central site 
FIG. 10. Dependence of the charge correlation function 
on the distance to the central site for a transverse exchange 
7 = 2. Correlations with electrons at gray (black) sites are 
shown in gray (black). Note the factor of ten difference in the 
vertical scale between the top and bottom graphs. Values for 
several ladder sizes are presented with different symbols as in 
previous figures. 

From the analysis of the correlation functions, we con- 
clude that there is a phase transition in the two-channel 
Kondo lattice when the transverse exchange parameter is 
larger than unity. Spin and charge correlation functions 
evolve from the ones typical of a conventional Kondo 
lattice to those of a d-wave Kondo ladder with gapless 
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spin and charge modes. The value of the spin correlation 
length derived from the scaling law (eq. (jfy) correlates 
only with the spin correlation length along the gray chain. 
The spin correlation length along the black chain appears 
to be determined by the hopping interaction rather than 
by the d-wave exchange coupling. On the other hand, 
the charge correlation lengths of the two zig-zag chains 
(black and gray) behave similarly and are in agreement 
with the lengths associated with the charge gap scaling. 
In any case, we have to keep in mind that results for 
larger ratios of Jd/Js are likely to be less reliable given 
that for these values 2£ ~ Ldiag, where Ldiag — 5 is the 
length for which an exact diagonalization is carried out. 



TABLE II. Charge-charge correlation lengths for different 
values of the ratio Jd/Js- The second column displays the 
lengths derived from the scaling law. The third and fourth 
columns show the correlation lengths between the central site 
and sites in the gray and black chains, respectively. 



can also be clearly seen in Fig. IO, where the solid line 
is the fitting function. 

In contrast with the case of the conventional Kondo 
ladder, the local and the conduction spins display very 
different behavior in the d-wave limit (see the bottom in- 
set in Fig. [ll]). The third component of the local spin of 
each rung is always positive; it decreases from the center 
((S§) ~ 0.44) to the end site of the ladder ((S£> - 0.16). 
On the other hand, the magnitude of the z-component of 
the electronic spin is always negative and ranges between 
—0.28 and —0.14. However, the sum of both spin densi- 
ties, (Sj total ) behaves similarly to the case of only on-site 
exchange coupling, although it is clear that the dispersion 
of the data is larger (Fig. [ll]). By using the same two- 
parameter fitting, {Sj total ) — C[sin 
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Let us conclude by discussing some characteristics of 
the spin excited states. Fig. [ll] displays the value of the 
third component of the total spin in each rung as a func- 
tion of the distance between the rung and the center of 
the ladder for Jd/J s = 0, open circles, and Jd/J s = 10, 
open diamonds. In the insets, the z-component of the 
spin of the local moments, {Sj), and the conduction elec- 
trons, (Sj), are displayed. The results correspond to a 
ladder of 12 sites with open boundary conditions. 

In the case of the conventional s-wave limit, the spin of 
the local moments, {Sj), and of the electrons, (s|), are 
almost identical for any site in the ladder, as it can be 
seen in the top inset of Fig. O. The magnitude of the 
spin is also equally distributed between the lower and up- 
per chains. The local and electronic spins follow-a square 
sine function, appropriate for a massive mode,c3 



{S 



j,total) ~ {Sj 
1 



J open 



sin 



") ~ 2<6y> 

n(j + (L 



2(4) 

open/2)) -i 



(7) 



J open 



where j is the position under consideration and L open is 
the length over which the spin excitation extends, which 
is longer than the ladder length. By making a two- 



parameter fit, {S* 



l ) = C[sm 



we infer 



J j, total/ — w L 0111 
-^open — 12.16 and a normalization coefficient C=0.0822, 
in agreement with the value of l/L open . This agreement 



7T(j + (L ope „/2)) 1 2 



we 



obtain a length of L open — 13.02 and a normalization 
coefficient C=0.0766, which should be compared with 
l/L open = 0.0768. 




2 3 4 

Distance to center of ladder 
FIG. 11. Third component of the total spin of each rung 
as a function of the distance to the center of the ladder 
for Jd/Js = 0, open circles, and Jd/Js = 10, open dia- 
monds. The lines represent the best fit to the expression: 
SLtai = ao(sin [ir(j + (L open /2))/ L open }) 2 . The solid line fits 
the Jd/Js ~ points, dashed line fits the Jd/Js = 10 set 
of points. The insets display the z-component of the spin of 
the local moments (triangles up) and the conduction electrons 
(triangles down) for these two values of the ratio between the 
exchange couplings. 



V. VARIATIONAL CALCULATION 

In order to gain more insight into the nature of the 
ground state in the d-wave limit, we have performed a 
variational calculation. Variational wave functions have 
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been successfully used in the conventional Kondo lattice 
modelxil. We restrict ourselves to a simpler Hamilto- 
nian with only d-wave coupling. In this case, as we have 
discussed previously, the ladder decoupled in two discon- 
nected chains (Fig. fij). We can label the sites of one 
of the zig-zag chains just by their j coordinate. So, the 
Hamiltonian for one of the decoupled chains is reduced 
to: 



2(2 



■££• 



(4,(-7K.'d, v (- 7 )) (8) 



with dj s (x) 



^ j+hs +ipj- hs + xip js , 



where ip^ s creates 



a conduction electron with spin s on site j. 

Our variational wave function is the product of a 
Zhang-Rice-like singlet at each site: 



l*o> = IIO TKxH-iD.d^^lo), 



(9) 



where | T)z (| J,);) indicates that the 1th spin is in the state 
with S z = +1/2 (S z = -1/2). The operator dj^a) cre- 
ates a non-orthogonal Wannier state which overlaps with 
the states at neighboring sites. We allow a to vary over 
values that differ from 7. 

By decomposing a into its magnitude and its phase, 
a = |a|e 1 ^, we find that the energy minimum for any 
transverse exchange parameter (7) and length of the 
chain is always reached when <f> = ir. Table [II displays 



the values of \a\ for which the variational energy becomes 
minimum and the value of that minimum. For 7 = 2, the 
asymptotic minimum is in reasonable agreement with the 
computed value of —1.05. The variational energy ob- 
tained for 7 = 0.5 is also remarkably close to the value 

obtained in the DMRG calculation (E™™° nd 0.88). 

However, the variational energy obtained for 7 = 1 is too 
high; it is even higher than the energy of the configu- 
ration when half of the spins form singlets with one of 
their nearest neighbors (E™ound = ~^/4). The energy 
minimum for 7 = 2 is a very shallow one, and there is a 
large region of \a\ values where the energy is very close 
to — 1. In contrast, for 7 = 1 and 0.5 the energy mini- 
mum is reached at finite values of the parameter \a\, and 
a scaling in 1/L is necessary to obtain the ground state 
energy of an infinite system. 



TABLE III. Results of our variational calculation for the 
studied transverse exchange parameters. The second and 
third columns show the magnitude of \a\ where the energy 
minimum is achieved and its value, respectively. 
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FIG. 12. Correlations of variational wave function ([]) 
versus the parameter |q|. Spin correlation between 
the central moment and the electron along the rung 
((S z oS z q)), solid line; correlation with the nearest elec- 
tron along the chain ({S^oS 2 ?}), dotted line, and corre- 
lation with the next-next-nearest electron along the rung 
({S z qS Z 2}), dashed line. The on-site charge-charge corre- 
lation, C(0) = (p(0)p(0)} - (p(0))<p(0)}, is also displayed 
(long-dashed line). 



We have also calculated the charge-charge and spin- 
spin correlations as functions of |a|, for <f> = ir. The 
spin-spin correlation between the local spin at the center 
of the ladder and the two nearest, and next-next-nearest 
electrons sitting in the black chain are displayed in Fig. 
|l2] , together with the on-site charge-charge correlation at 
the central spot. Since gray and black chains are totally 
decoupled in our model (Hamiltonian (^|)), correlations 
with electrons at the gray sites are absent. From this 
plot we conclude that extended Kondo models, where the 
moments couple with electrons on other (rather than it's 
own) sites, display extended spin correlations and strong 
on-site charge-charge correlations. There is a large range 
of values of the variational parameter \a\ where the spin 
and charge correlations reach big values. For example, 
for \a\ = 2.8, C(0) = 0.24 and (S z0 s zl ) = -0.18. These 
values are in good agreement with those found in our 
DMRG calculation for 7 = 2. So, it is reasonable to 
think that the ground state of the d-wave Kondo ladder, 
at least for some values of the transverse coupling, can be 
understood in terms of extended singlets on each lattice 
site (Eq. |). 







VI. CONCLUSIONS 

Motivated by experimentally and theoretical consid- 
erations—About the phase diagram of Kondo lattice 
modelsJEI we haver-used the density matrix renormaliza- 
tion group methodE3 to study the strong coupling limit of 
a Kondo lattice model with two different interactions: the 
conventional on-site Kondo exchange (J s ) and a d-wave 
exchange coupling (Jd)- Our purpose has been to map 
the phase diagram of this model and sort out whether or 
not it displays a phase transition. 

By increasing the ratio between the two couplings 
Jd/Js-, the system smoothly evolves from a conventional 
Kondo insulator ground state, where a local singlet is 
formed at each site, to a ground state with much smaller 
spin and charge excitation gaps. For small values of 7 less 
than unity, our results suggest that the system preserves 
a spin and charge gap for all values of the Jd/ J s ratio. 
However, when the parameter 7 is increased to values 
larger than unity, there is clear evidence to suggest that 
the system develops gapless excitations for Jd/J s ~ 10, 
suggesting a phase transition. For a transverse exchange 
parameter of 7 = 2 and within the accuracy of the cal- 
culations, both the spin and charge gap become zero at 
this point. However, numerical difficulties related to the 
proximity to an antiferromagnetic instability preclude us 
from resolving definitively whether these gaps are zero 
for Jd > 10J S corresponding to a phase transition into 
a metallic phase, or whether the gaps re-establish small 
but finite values at Jd > 10 J s , corresponding to a critical 
point at Jd/Js = 10. Our intention is to increase the 
Hilbert space and the length of the ladder to be able to 
get more accurate values of the excitations gaps in this 
limit. 

This semimetallic ground state with gapless spin and 
charge excitations can be understood as a superposition 
of the extended singlets formed between each local mo- 
ment and its three nearest-neighbor electrons. The char- 
acteristics of the ground state excitations change with 
Jd/Js- In particular, we find that as Jd/Js rises, the ra- 
tio of charge to spin gap A c /A s is enhanced. Simultane- 
ously, the system develops strong spin-spin correlations 
between the local moment and the electron along the 
rung and large on-site charge-charge correlations. The 
typical correlation lengths increase with the Jd/Js ratio. 
Also the transition from an exponential to a power law 
decay in the spin and charge correlations for 7 = 2 and 
Jd/Js > 10 supports the existence of a gapless phase. 

Additional insight into the results at large Jd has been 
obtained from a variational calculation. By writing down 
a variational wavefunction in which an extended singlet 
forms between each moment and the three neighboring 
electrons in the ladder, we obtain variational ground 
state energies that are close to our DMRG values for 
7 = 2 and 0.5. Our variational ansatz also displays large 
spin-spin correlations along the rungs and strong on-site 
charge-charge correlations as found in the DMRG results. 



It would be interesting in future work to examine the 
effect of three-leg and four-leg ladders on the above re- 
sults. By increasing the number of legs, the full point- 
group symmetry of the two-dimensional limit is gradually 
restored, permiting d-wave singlets to delocalize without 
admixing into extended s-wave pairs. This may permit 
the gapless phase to be reached at lower values of the 
ratio Jd/Js- Indeed, the vanishing spin and charge gaps 
observed in our ladder model may become a full-fkdged 
phase transition in truly two-dimensional systems.Eil 
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